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CN . ABSTRACT 

' The magneto-rotational instability originates from the elastic coupling of fluid elements 

Q , in orbit around a gravitational well. Since inertial accelerations play a fundamental 

dynamical role in the process, one may expect substantial modifications by strong 
gravity in the case of accretion on to a black hole. In this paper, we develop a fully 
, covariant, Lagrangian displacement vector field formalism with the aim of addressing 

these issues for a disk embedded in a stationary geometry with negligible radial flow. 
This construction enables a transparent connection between particle dynamics and 
the ensuing dispersion relation for MHD wave modes. The MRI-in its incompressible 
variant-is found to operate virtually unabated down to the marginally stable orbit; the 
\ putative inner boundary of standard accretion disk theory To get a qualitative feel for 

■ the dynamical evolution of the flow below r ms , we assume a mildly advective accretion 

' flow such that the angular velocity profile departs slowly from circular geodesic flow. 

£NJ , This exercise suggests the turbulent eddies will occur at spatial scales approaching 

' the radial distance while tracking the surfaces of null angular velocity gradients. The 

implied field topology, namely large-scale horizontal field domains, should yield strong 
mass segregation at the displacement nodes of the non-linear modes when radiation 
I 1 stress dominates the local disk structure (an expectation supported by quasi-linear 

arguments and by the non- linear behavior of the MM in a non-relativistic setting). 
Under this circumstance, baryon-poor flux in horizontal field domains will be subject 
to radial buoyancy and to the Parker instability, thereby promoting the growth of 
poloidal field. 
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1 INTRODUCTION 

The process of accretion onto compact objects has been long recognized as the primary mechanism to power the most luminous 
events in space. In the traditional picture of Shakura & Sunyaev (1973) and Novikov & Thorne (1973), entropy is generated 
and radiated locally from the free energy available in a shear flow with a Keplerian angular velocity profile. Two salient 
oversimplifications of this framework have been the focus of intense research and progress in the last decade: energy advection 
by the flow and free-energy tapping and angular momentum transport through magnetohydrodynamical processes. 

The magnetorotational instability or MRI (Velikov 1959, Chandrasekhar 1961, Balbus & Hawley 1991), justifies the 
long-sought mechanism for efficient, turbulent transport of angular momentum that enables accretion disks to operate with 
astrophysically interesting mass accretion rates (Pringle 1981). The importance of this process cannot be overstated: By 
catalyzing accretion into gravitational wells, the MRI enables a plethora of astrophysical phenomena to occur, from protostar 
formation inside molecular clouds to jet launching in quasars. The MRI also holds the key to understand the extraction of 
free energy from the differential shear flow of otherwise hydrodynamicalty stable disks (Balbus et al. 1999, Godon & Livio 
2000). 

On the observational front, the wealth of high-quality data from spectral and timing devices aboard space-borne high- 
energy observatories has turned out the most compelling evidence yet of accretion onto black holes. The discoveries of pairs of 
high frequency quasi-periodic oscillations in RXTE X-ray timing data from microquasars GRO J1655-40 and GRO 1915+105 
(Strohmayer 2001a,b) have brought the spotlight to hydrodynamical models of adiabatic global excitations of the inner disk, 
a.k.a. diskoseismology models (Perez et al. 1997) or relativistic precession models (Stella et al. 1999). Interestingly, the role 
of magnetic fields has been largely ignored in spite of clear evidence that QPO's, being non-thermal, hard X-ray phenomena, 
likely do not originate in the accretion disk proper but rather on a magnetically active accretion disk corona (R. Blandford, 
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Priv.Comm.). Likewise, the recent report (Wilms et al. 2001) of the detection of a very broad Fe Ktv feature on the XMM- 
EPIC spectrum of MCG-6-30-15, has made a very strong case for the inadequacy of standard models of energy deposition 
in accretion disks. The proposed exits to this paradox-extraction of black hole spin energy (Blandford & Znajek 1977), or 
non-zero torque at the marginally stable circular orbit radius, r ms (Agol & Krolik 2000)- both rely on magnetic coupling 
between a standard disk and the flow inside r ms . 

On a more exotic front, theoretical progress in our understanding of accretion processes at the most extreme imaginable 
conditions-stellar-mass black holes hyper-accreting at twelve orders of magnitude above the Eddington limit-requires attention 
to a detailed physical account of highly relativistic accretion flows. Aside from the potential to explain gamma-ray burst 
phenomenology, such studies are chiefly relevant to assess the likelihood of "failed" supernovae as gravitational wave sources 
(Fryer et al. 2001). Indeed, when neutrino trapping occurs at M^IMq sec -1 (Popham et al. 1998), the associated dynamical 
stress will mimic the effects of radiation stress in standard disks where clumpy accretion ensues (Turner et al. 2001). If the 
mass fraction in the clumps is large, prolific gravitational waves will be emitted from the mass quadrupole moment associated 
with the bulk motion of large mass over-densities. Such a scenario will also lead to excitations of the black hole's geometry 
which, at high values of the spin parameter a, can produce highly characteristic, monochromatic black hole ringing as the 
geometry settles towards a quiescent Kerr state (Araya-Gochez 2003). Remarkably, the expectation of a large mass fraction 
in the clumps is reasonable and justifiable by the physical picture of near-hole accretion presented herein. 

An outstanding issue yet to be addressed in light of recent theoretical progress is our view of black hole accretion inside the 
marginally stable orbit, the putative inner boundary of standard accretion disk theory. In particular, very little is concretely 
known about the inertial effects of strong gravity on the relevant MHD processes. Previous work has either assumed pure 
hydrodynamical flow (and energetically negligible energy release) or, alternatively, laminar flow under ideal MHD conditions 
(Krolik 1999, Gammie 1999). Krolik (1999) has brought out an interesting point: Under mere flux freezing conditions the 
assumption of ballistic orbits in the plunging region is never self-consistent; when the radial velocity component is significant, 
the magnetic field energy density becomes comparable to the rest mass energy density of the matter. 

In this paper, we address the issue of stability of the magnetic field (co-moving frame) in a stationary, axially symmetric 
background geometry. Curiously, the two key developments in accretion disk theory of the last few years may have come of age 
to properly address the problem at hand: Inside r ms the accretion flow will be mildly advective, with a slightly sub-Keplerian 
angular velocity profile and possibly supported in part by the radial pressure gradient of a hot MHD fluid with significant 
relativistic enthalpy (see Popham & Gammie (1998) solutions for moderate values of a and advected fraction /). In this spirit, 
we argue in §5 that the natural evolution of the MRI inside the marginally stable orbit is at least consistent with this view. 

The (magnetohydro) dynamics of black hole accretion comprises two important aspects that have received relatively 
little attention: the effects of radiation pressure (see, however, Blaes & Socrates and Turner et al. 2001), and the effects of 
strong gravity (see footnote 7 of Gammie & Popham 1998). We will address the former problem in a future paper (Araya- 
Gochez & Vishniac 2002) while concentrating on general relativity in this one. As a background, §2 looks at the Lagrangian 
displacement vector field formulation of the MRI concentrating on inertial and compressibility effects. In section 3, we develop 
a fully covariant theory of the instability. The intention is to build a theoretical framework from first principles in order to 
avoid missing any subtleties associated to the full incorporation of gravitational effects (e.g., reference is made to the Cowling 
approximation and to the fixing of the gauge associated with the component of the Lagrangian diplacement along the fluid's 
four- velocity). The elastic reponse of the field is computed by noting that the surface of invariance of the Faraday tensor 
attributes mathematically identical variational properties to the two four- vectors that span it: the magnetic field four- vector 
and the fluid's four-velocity. We then make the minimal modifications to the relativistic fluid equations that allow for the 
inclusion of a coherent magnetic field and undertake a local stability analysis of this field in the medium of a slim disk around 
a rotating black hole, while suppressing compression. The role of compressibility in a photon gas is then briefly assessed. 



2 A LAGRANGIAN FORMULATION OF THE MRI IN COMPRESSIBLE MEDIA 

The MRI is essentially a local instability. In the frame of the fluid, the interplay of inertial "forces" with the elastic coupling 
of fluid elements creates an unstable situation to the redistribution of specific angular momentum, I. Without the elastic 
coupling provided by the bending of field lines, such inertial forces-namely, the shear(tide) and the coriolis terms-induce 
radial epicyclic motions while preserving specific angular momentum in collisionless fluids (e.g., stars in the Galaxy). This is 
related to the Rayleigh criterium for stability of a differentially rotating fluid: r~ 3 d r £ 2 — k 2 > 0, where n is the frequency of 
epicyclic motions. 

In the weak field limit, one may construct a dispersion relation quite independently of the specific magnetic field topology: 
Highly sub-thermal fields, v AU /c s <JC 1, guarantee that the instability is truly local 1 , occurring at large values of fey ~ Q/v AU (= 
k- ljj). In this simplified approach, the global disk structure is ignored (no curvature nor radial structure) and the response of 
the field amounts to nothing more than providing a restoring force to displacements from equilibrium (Balbus & Hawley 1992). 
Indeed, in the horizontal regime of Lagrangian displacement two orthogonal field topologies yield nearly identical mathematical 
dispersion relations for wavemodes: axisymmetric perturbations of a meridional field and non-axisymmetric perturbations of 

1 When the field is non-negligible, fc^ 1 may approach the disk's pressure scale height and in the case of supra-thermal toroidal fields, 
non-axisymmetric modes have fastest growing wave numbers that may approach the inverse radial scale- length (Foglizzo & Tagger 1995). 
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a toroidal field. The former case corresponds to the "classical" Balbus-Hawley instability and it's physical relevance is free of 
controversy. The relative importance of the latter analyses is a more subtle issue. 

A somewhat technical point-well discussed in the review by Balbus & Hawley (1998)-is the non-locality induced by 
shear on wave-modes with Eulerian coordinate phase-dependencies. For k v ^ modes, shear evolves the radial component of 
wavenumbers according to k r (t) — ko r — [di nr Q] k v t, which means that modes that could be "unstable" are only so, transiently. 
The maximum instantaneous growth rate occurs when k r — > and matches that of the local axisymmetric modes. In the 
end, this issue turns out to be more academic than practical but it stresses the importance of treating the instability locally, 
in co-moving coordinates. The down side is that this greatly complicates global approaches that rely on eigenmode solutions 
in Eulerian coordinates extrinsic to the fluid. On the other hand, in a local approach azimuthal wavenumbers are no longer 
discrete (Ogilvie & Pringle 1996) and consequently, neither are the co-moving frequencies (see below). 

A related issue concerns the relevance of non-axisymmetric mode analyses when the magnetic field is not purely toroidal. 
Balbus & Hawley (1998) argue that the strict ordering of wavenumber components (and narrow phase space) necessitated to 
achieve fastest growth: k v <ti k r <JC ke , ensue in violent poloidal Alfvenic couplings that promptly take over the dynamics. Non- 
axisymmetric modes, however, are important for at least two key reasons: a- the ordering is not so restrictive when the fields 
are not weak (as needed to explain a values of a few tenths), and b- compressive, non-axisymmetric modes are fundamental 
to examine energy deposition when radiation stress becomes significant (Araya-Gochez &: Vishniac 2002). Moreover, because 
the dispersion relations relate simply (at least in the horizontal regime of fastest growth), it is rather lucrative to examine 
both cases at once. 

Aiming to formulate a fully covariant relativistic theory of the MRI in §3, this section conducts the same in three 
dimensions. The linear stability analysis is carried out in terms of the Lagrangian displacement vector field, £. Foglizzo (1995) 
has stressed the usefulness of this approach to account for the polarization of compressive MHD modes. A simple meridional 
stratification profile sets the physical scale-length of the problem: d z hip = HT 1 , with gas, radiation (and possibly magnetic) 
pressures tracking the unperturbed density profile pHSl = jJ r+g +P B - The problem naturally splits in two parts: computation 
of the inertial/geometric terms §2.1, and computation of the body forces from gas, radiation and electromagnetic stresses 
§2.2. We avoid going into the rotating frame from the onset in order to preserve a transparent connection to a "universal" 
standard of rest frame (to be associated with Boyer-Lindquist coordinates). 



2.1 inertial terms 

Inertial accelerations are geometrically imprint in the connection terms for the covariant derivatives of the Eulerian velocity 
components. For spherical coordinate motion (r,ip,6) — > (V r ,V V ,V 6 ), the only non-trivial connections are V r vv A Vf v . 
Denoting the Lagrangian time derivative by d t = d t + V • V, the three components of Euler's equation read 

d t V r = d t V r + V J 'Vj + {-r)V p V' p = g rj fj 

d t V v = d t V v + V j V? + {2/r)V r V v = g^fj 

and d t V e = d t V e + V j V% = g ej ij 
1 1 

where f=--Vp+- — J x B, (1) 

p A-Kp 

g %i is the flat-space metric for spherical coordinates, and V will stand for the covariant derivative hereon. 

Assuming an equilibrium from purely azimuthal (but differential) bulk motion V= an Eulerian perturbation of such 

a state (V r , V v , V z ) — > (v r , Q. + v v ,v z ), leads to the usual equations associated with a rotating frame and its coriolis and 
centrifugal terms. For coordinate motion, Euler equations for the perturbations of the fluid read 

{dt + Vtd v )v r - 2rtlv v = g rj Sij 

{d t + 0.d v )v e = g ej 5fj (2) 

where 8f stands for the Eulerian perturbation of the sum of specific body forces. The standard form of these equations, e.g., 
for non-coordinate motion (see, Chandrasekhar 1961), may be obtained from Eqs [1] above by "dimensionalizing" V v , (i.e., 
in the second Eq, multiplying by r and completing the differential while recalling that the covariant derivative and the metric 
commute [V, gij] = 0). 

Next, one switches dynamical variables from the Euler velocity perturbation, v, to the Lagrangian displacement, £, using 
the first order relation between Lagrangian and Eulerian variations, A = 5 + £ ■ V, whilst denoting 2 AV = dt£ and <5V = v 
(see, e.g., Chandrasekhar & Lebovitz 1964, Lynden-Bell & Ostriker 1967) 

2 The tilde indicates that this form of Lagrangian displacement-which is generally non- unique-has had its gauge "fixed" in accordance 
to the non-relativistic regime. Mathematically, this amounts to a choice of Universal time direction, It (e.g., unaffected by the fluid's 
motion), while adopting the gauge fixing condition £ • It = 0. 
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v = {d t + v ■ v} £ - (| ■ v)v — icr* - f n, r i„. (3) 

The algebraic relation follows from the assumption of differential rotation and from writing exp i(u>t + m<p + k z z) dependencies 
for Note that the connection coefficients in Eq [3] cancel one another and that a = co + mQ denotes the co-moving frequency 
of the perturbations. 

These geometrical equations have their more traditional equivalents in the so-called shearing sheet approximation where 
a co-moving, "locally Cartesian frame" (r,ip,6) — > (x,y,z), is used along with the linearized shear velocity field, V(a;) = 
[di nr f2] xl y , to treat the problem locally while introducing the coriolis terms by hand. Defining the Cartesian derivative 
operator d, then the equivalent to Eq [3] is (dt + V • d)£ = v + £ • 9 V, which has Galilean invariance in the sense that 
Lagrangian time derivatives produce co-moving frequency factors in the dispersion relation: (d t + V • d)£ = c? t £ — > ia£. 

The equations of motion (EoM) for the (coordinate) Lagrangian displacement are 

(-a 2 + 20. rSl, r )Z T - 2rVL ia £ v = Sf r 

-a 2 ^ + 2-iaC = Sf v 
r 

-a 2 ^ = 5f. (4) 

Note that the Eulerian shear term, oc l v , becomes the tide term, oc l r , in terms of £. 
Let us re-cast these equations in a more compact form 

e + 2r jk v j £ k - 2v) k v ] {v - i) k = sr (5) 

where each over-dot stands for a factor of ia (from a Lagrangian time derivative). In the shearing sheet approximation, 
these equations correspond to the Hill equations for non-coordinate motion (Chandrasekhar 1961, Balbus & Hawley 1992) 
£ + 2 f2 x £ + 2rQn, r ^ l x = 5f . 

2.2 compressibility 

The Lagrangian perturbation of mass density and the Eulerian perturbation of the field follow from mass and magnetic flux 
conservations (recall the non-relativistic relation A = 8 + £ ■ V) 

^ = -V-£, <5B = Vx(£xB). (6) 
P 

The latter equation includes possible gradients of the background field 9B =fc 0; however, in the spirit of examining the 
instability as a local phenomenon the global structure of the field is ignored herein. The Lorentz force variation (co-moving 
frame) may then be written as 



^Jxb) = V f x { V(V.0+V|£-V 6 V(1 6 .£)-1bVb(V-0} 



~<u x - fc B £ B )k + 4€ - (7) 

where the scalar operator Vg = 1 B • V, and lg is a unit vector in the direction of the unperturbed field. Note that the term 
(ki£,i — fc B £ B ) = k±£± may be interpreted as a restoring force due to the compression of field lines (distinct from line bending, 
Foglizzo & Tagger 1995). 

The Lagrangian variation of the specific pressure gradient contains two terms (Lynden-Bell & Ostriker 1967): one oc Ap _1 
and another oc AVp r+g - In terms of the displacement vector, the first term is proportional to the equilibrium value of Vp r+g 
which is negligible 3 in the local treatment (proportional to a radial gradient). For the same reason, the Eulerian and Lagrangian 
variations of the pressure "force" are identical. 

The thermodynamic pressure term is then given by 

-*(iv Pr+g ) =r^±iv(v.£) v ^ c 2 kk^, (8) 

where, for heterogeneous media, F = d[i np ] lnp represents a generalized adiabatic index (see, e.g., Chandrasekhar 1939, Mihalas 
& Mihalas 1984). 

Putting the above equations together, one gets the EoM in Fourier-space 

e + 2r} fc V^ fc - 2T) k Vi(v - i) k = g^[{c 2 s k + (k - fc B l B )}(k ■ |) + < f (4 | - fc B C B k)fc. (9) 

which agrees with the matrix de-composition of Foglizzo & Tagger (1995) in the case of a purely toroidal field embedded in 
a gas with adiabatic index F — 1 . 

3 the variation of the mass density is also negligible when the focus is on the effects of radiative heat conduction: Loss of pressure support 
out of compressive modes involves only the pressure term (Araya-Gochez & Vishniac 2002). 
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Reckoning of fluid compressibility has complicated somewhat the equations of motion. Yet, these generally unwieldy 
equations simplify greatly in the regime of fastest growth (a.k.a. the horizontal regime) and for two ideal field topologies of 
interest. When the field is meridional, the fastest growth modes have £g = 0, and k ~ fcglg thus yielding a simple isotropic 
elastic response oc — v^ if (ik^) 2 

Alternatively, when the field is purely toroidal, the meridional component of Eq [9] yields an anisotropy constraint: 
v ah (^x£±) ~ — c ? (k ' £) (Foglizzo & Tagger 1995), which allows for a straightforward solution in this regime. Defining A 
through 

1 " A = "5| = f^0 < 10 > 
where = p Bv /p r+g , the dispersion relation out of Eq [9] reads 

o A - { (A + l)g| + X 2 } a 2 + A q\ {q\ + 4A} = 0, (11) 

where all frequencies are normalized to the rotation rate, A = idi nr hif2 is the Oort A "constant", \ 2 = 4(1 + A) is the 
squared of the epicyclic frequency, and q & = (k • VAif ) /fi is a frequency related to the component of the wave vector along 
the field (in velocity units). 

The non-axisymmetric modes of fastest growth conform with (Araya-Gochez & Vishniac 2002) 



2A + (i±A)x {-^-}, where V = 1 + (i=±)A + ^ 1 + (1 - A)A, (12) 



! 2 

V 

while the expression in the curly brackets identifies with the negative root of the dispersion relation. 

Note that the compressibility of non-axisymmetric modes is imprint on the deviations of A from unity. From Eq [10] one 
reads that the degree of compression of these modes gets stronger with the (toroidal) field strength and, naturally, with a 
softer equation of state. In the companion paper (Araya-Gochez & Vishniac 2002), we find that when radiation pressure begins 
to dominate the disk dynamics, an "ultra soft" effective (adiabatic) index accentuates the effects of mode compressibility. On 
the other hand, setting A = 1 and re-orienting the field vertically produces the standard (incompressible) dispersion form for 
the Balbus-Hawley instability of a meridional field in the horizontal regime. 



3 GENERAL RELATIVISTIC EFFECTS IN THE COWLING LIMIT 

In contrast to the Newtonian case, the formulation of a covariant theory of accretion disk oscillations requires more than 
mere application of the Lagrangian rate of change operator d t (or its relativistic counterpart d T ) to the Eulerian velocity 
perturbation. This is insufficient to carry out a normal mode analysis because of the freedom associated with the choice of 
coordinates. It is much more useful and proper to free the eigenmodes from the coordinate representation; treating them 
rather as intrinsic to the physical system. It is here that a Lagrangian construction comes in handy. 

A working covariant definition of the Lagrangian displacement is that of a vector field that moves a fluid element's 
world line from its unperturbed position in spacetime to its perturbed one. The fundamental relation between the Lagrangian 
(following the fluid's world line) and Eulerian (taken at a fixed coordinate point) variational operators is 

A = 5 + C^ (13) 

where £ stands for the Lie derivative. 

An elemental use of this relation involves particle number conservation (Schutz & Sorkin 1977): with the use of a number 
flux density A/"" = n v / — g U" , such a law reads AM = 0, where g = det \g^ v \ and U" is the four- velocity of the fluid. In the 
Cowling approximation, Sg = 0, and in the absence of co-moving sources (or sinks) of particles, one gets for the variations of 
the four-velocity of an ideal fluid: 

AU = = SU + C ( U (14) 

which demonstrates that Eulerian perturbations of the four velocity, SU = u, obey u v — —C(U V . 

The connection to the Newtonian limit is recuperated upon identifying c U ■ V with the convective (or material) rate of 
change cU ■ V c — ^? (dt + V • V) so that, with AV = d t £, one has 

A = 5 + £-V (15) 

(see, e.g., Chandrasekhar 1964, Lynden-Bell & Ostriker 1967, and compare SV = v with Eq [3]). 

The fluid particles that constitute a thin accretion disk (with negligible radial inflow) embedded in a Kerr spacetime 
geometry have unperturbed four velocity U v = 7(1* + where l t = (1, 0, 0, 0) and \ v = (0, 0, 1, 0) are the Killing vector 

fields of the stationary, axisymmetric geometry and where 7 is the "redshift" factor of the fluid elements at fixed radius 
7 = U = d T t. In terms of the Lagrangian displacement vector field, each Lie derivative with respect to one of the Killing 
vector fields of the geometry "brings down" a wavenumber co-factor in the dispersion relation (modulo spatial gradients of 
the four velocity) . This leads to an algebraic relation between £ and u = 5U in the case of a differentially rotating fluid: 

u = £ 7( i t +«i„)£ 
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= £-7£ r n, r l»,-C^-Vln7. (16) 

Here £ = iayt^ , with a = + mf2 the co-moving frequency of the perturbation as measured at asymptotic infinity (see Ipser 
& Lindblom 1992). 

The relativistic generalization of Eq [3] is found upon projecting £ on the 3-surface perpendicular to the (unperturbed) 
four velocity. With h al3 = V^U 13 + g al3 the projection operator, one has 

h a p u fi = u a = ia-yC -yf Q, r i° (17) 

while fixing the gauge freedom associated with the component of the Lagrangian displacement perpendicular to space-like 
hypersurfaces (Schutz & Sorkin 1977), i.e. along the local "time" direction. Note that the requirement of unit normal for the 
perturbed velocity, U + u, under the Cowling approximation fixes the gauge accordingly: 2U a u a = — U a U^ Sg a/ 3 = 0. 

Dynamical conservation laws for an ideal fluid in the presence of a large-scale electromagnetic field are written succinctly 
through the Einstein-Maxwell equation (Cowling approximation) 

T^" v — F^Jv = 0, 

where the first term stands for the matter stress and the second equals the Maxwell stress. The notation is standard fare: F is 
the electromagnetic field tensor and J — neU is the four-current. Ideal MHD makes things easy by stating that the electric 
field in the co-moving frame vanishes everywhere. Since the latter is the contraction of the field tensor with the four-velocity, 
it follows that F ilv .l v = 0, and the four-acceleration from the Maxwell stress vanishes as well (but not its perturbation). 
We shall concern ourselves with the material stress first. 

Denoting the relativistic enthalpy of the fluid by g = p + e + p, it is straightforward to show that for a non-dissipative, 
ideal fluid such that T M " = gU^U v + pg^ , 

Tr„ = gA T U" + g fl "p,„ where d T = U • V (18) 

stands for the generalization of the convective rate of change, i.e. the Lagrangian proper-time derivative. The projection of 
this equation along the four-velocity states energy conservation while the perpendicular components express conservation of 
momentum. 

The specific Eulerian perturbation of the four-acceleration (normalized to the enthalpy) looks like «-V U + L^V u, and one 
can use Eq [17] to switch the dynamical variable in favor of the projected £'s 4 : 

u-X7U + h(U-X7u a ) = h(U-V£-£-VU)+2£-VU 

- 7 £ r firVVL/-ft[/-V(7r«ri ¥ >) (19) 

where the silly hats on the £'s (signifying projected components) were dropped. 

The first term on the r.h.s. of Eq [19] may be readily identified with (the projection of) the Lie derivative of £ along U. 
Defining q = l t + (so that U = jq), one computes 

Cui = 7A£-?£-v 7 

= £- 7 £ r fVi ¥ , + £t/-Vln7-£/fVln7 (20) 

where £ = (iuj) 2 ^. Note that the last term disappears upon (re)projection onto proper space-like hypersurfaces. 

The second term on the r.h.s. of Eq [19] is easily evaluated, £-V U a = j^Q^l^ + r™„{y M £", and the third simply involves 
a projected affine connection. Evaluation of the (non-projected) last term yields four pieces: 

[/■V(7Cfiri v ) = C^,rhU-\7-y + 'yQ,ri v U-\7C 

+ 7 £% U-VSl,r + 7C",rt/-VV 

Under the premise of negligible radial motion, in the second piece above C/-V£ r ~ £ r + 0(U r ), while the third is 0(U r ). 
Likewise, pairing of all terms proportional to the logarithmic gradient of the redshift factor yield same order (negligible) 
corrections (£ — y^Q^l^U-V In 7 ~ 0(U r ). Moreover, the last term above involves the same connection coefficient as the 
third term on the r.h.s. of Eq [19]. 

When all this is said and done, one gets for the specific Eulerian perturbation of the four-acceleration: 

u-vu + h (u-v u) ^-^ t + 2r^tr^ - 2r v afj u a {u - if. (21) 

The resemblance with Eq [5] is remarkable but not accidental. 

The shear(tidal) term is embodied by the third term on the r.h.s.: u — £ = — 7 £ r f2 jr l v . Note the non-trivially hatted unit 

4 note that with this form of the stress-energy tensor, h a ^d T U^ = d T U a , i.e. the four-acceleration automatically lies in proper space-like 
hypersurfaces. 
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vector = h^V^p = 1£, + U V U V . We evaluate this term first using the standard form of the Kerr metric in the equatorial 
plane (Boyer-Lindquist coordinates): 

ds 2 = ~^dt 2 + r 2 A(dip - u>dt) 2 + ^dr 2 , (22) 

with lo = 2a/ Ar 3 the rate of frame dragging by the hole and where the metric functions of the radial BLF coordinate are 
written as relativistic corrections (e.g. Novikov and Thorne 1973): 

A = 1 + a 2 /r 2 + 2a 2 /r 3 , and V = l-2/r + a 2 /r 2 , 

in normalized geometrical units (c = G — M bh =1). 

In expanded form, the projection of the Killing vector associated with the azimuthal symmetry is 

i v = [l + j 2 fv^n]l^+j 2 fv^l t , 

where 7 = 'yW'D/A is the redshift factor relative to "locally non-rotating observers" (Bardeen, Press & Teukolsky 1972) and 
f = rA/VV is the radius of gyration for the physical velocity in that frame, = f(Q — oj). 
Evaluation of the tidal term is a bit lengthy but straightforward 



-2^1^(6-0" = ~ {^ 2 vd lnr n} 



(23) 



where fi± = ±(r 3 ^ 2 ± a) -1 refer to prograde and retrograde circular orbits and where the expression in the curly brackets 
equals (minus) the shear of the congruence of circular, equatorial geodesies (Novikov and Thorne 1973). 

Next, to evaluate the coriolis terms one finds from the gauge fixing condition £ • U = 0. Accordingly one finds 

OT ,r 1jai D V ( » r 2 A(n-io) \ - 



and 



2r t ag U a £> 3 = 2~,-± 5 (r 2 + a 2 -an(3r 2 + a 2 )) t. (24) 

Eq [21] for the (Eulerian) perturbation of the four-acceleration, a M = w-V + h (f7-V -u 1 "), was derived for the components 
of the Lagrangian displacement in a coordinate frame that is fixed with respect to distant stars, i.e. in the Boyer-Lindquist 
"frame". However, because the instability is local (at least for weak fields in thin disks) one needs to transform the compo- 
nents of Eq [21] for manipulation in terms of the local tetrad carried by co-moving observers. This simply involves (matrix) 
multiplication by the basis vectors of such a tetrad (e.g., Novikov & Thorne 1973). In our notation, the relevant basis vec- 
tors are e r a — l/\/V(0, 1, 0, 0) and e% = 7rV2?(— fi, 0, 0, 1) (note that transformation to the local tetrad yields equations for 
non-coordinate motion, i.e. equivalent to motion in a local Cartesian basis). 

Transformation of the r-component is trivial (since one needs to transform both the acceleration and the displacement 
vector in the basic EoM below, the radial scale, 1/vP, has no net effect): 



r 

4 



{±j 2 vd lnr n} 



f- (25) 



On the other hand, using the local azimuthal base vector, Eqs [24] and the gauge fixing condition £ • U = 0, computation of 
the local (^-component, oc — Qa l + df , is a bit more involved (again, the radial scale factors out in the EoM and does not 
affect the dispersion relation) 

1 & (\ - 2(1 - 2aQ) - r 2 AQ, 2 \ - w 1 / fi , N , , 2s 



Let us take a look at the Maxwell stress next. 

The fundamental premise of ideal MHD may be stated rather succinctly: in the rest frame of the fluid currents will flow 
uninhibited to (instantaneously) cancel any hint of an electric field. A relativistic generalization of ideal MHD may be achieved 
by a similar covariant (albeit imperfect) postulate: E r .f. = F • U = 0, e.g., the Faraday field tensor is "purely magnetic" in 
the fluid frame. 

Such postulate brings a few mathematical consequences (Phinney 1983): 

i- the second electromagnetic invariant vanishes everywhere, 

\F^ V = (27) 

ii- the Faraday field tensor is Lie transported along the worldlines of the fluid, 
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C V F = (28) 

Hi- the other zero eigenvector of the field tensor is the (space-like) magnetic field B = T ■ U with T^" = e MI/a/3 F a/ 3 the dual 
to the field tensor, 

F-{JF-U) = ^ (29) 
and iv- the field tensor is also invariant when transported along the magnetic field four-vector: 

<^ T . V F = % (30) 

Note further that, since the four-velocity and the four-magnetic field are orthogonal, U ■ B = 0, properties ii and iv above 
define a 2-surface of invariance for the Faraday tensor 

£ a u +bB F = ® (31) 

where a and b are arbitrary real numbers. 

Aside from the intrinsic (physical) difference in their spacetime orientation, the mathematical similarities between U and 
B are uncanny. 

Let us go back to the Newtonian case for a moment. From our definition of the Lagrangian variation of the three-velocity: 
d t £ = AV, one finds the equation governing the Lagrangian change of three- velocity: AV = (V-V) £, (e.g. Eq [15]). As noted 
in the footnote, the difference between A and A is related to the choice of gauge for £ Q . The induction equation of non- 
relativistic MHD yields a virtually identical relation for the magnetic field variation (in a frame where the fluid was originally 
at rest) which is spoiled by fluid compressibility: AB = (B-V) £ B V-£. Nevertheless, making use of the continuity equation 
and weighting the field by the inverse of the fluid's mass density B = B/p cleans up its connection to the displacement vector 
field AB = BV£. If only conservative forces, U • f = 0, act on the fluid, it can be shown that use of energy conservation in 
lieu of mass conservation simply swaps the fluid's rest mass density by the relativistic enthalpy (a world scalar), above. Thus, 
we choose to work below with a specific measure of the magnetic four-vector weighted by the inverse of the fluid's relativistic 
enthalpy B = -{JF-U). Such combination of observables (and its perturbation) occurs naturally in the problem at hand. 

Applying the Lagrangian variational operator, c.f. Eq [13], on B under the constraints from ideal MHD noted above, Eq 
[31], yields (contrast this with Eq[14]) 

AB = = SB + C t B. (32) 

This equation states a manifestly covariant expression for the Eulerian perturbation of the (enthalpy-weighted) Faraday tensor 
under ideal MHD constraints: 6 = Demanding that the total magnetic field four- vector be orthogonal to the (unperturbed) 
four- velocity is equivalent to projecting its Eulerian perturbation into proper spacelike hyp ersurf aces: 6 M — * b M = h^C B ^. 
Again, we suppress the hats below while tacitly imposing the condition £ • U = throughout. 

The Eulerian perturbation of the specific measure of the Lorentz force, 5(T e m ) M £ = SF^ V J v + F^ V S3 V , may now be 
written in terms of the Lagrangian displacement but the general expressions are not particularly illuminating. Evaluated in 
the frame where the fluid was originally at rest, the Eulerian perturbations of the field tensor and of the four-current depend 
linearly on the components of b: 8F = F(b) and 4ir8J = d ■ SF(b). 

We proceed by assuming negligible gradients of the background specific field (VB = 0): 

fo = B-VC-C-VB — >B-(ik)£, (33) 

and, consistent with this assumption, we also ignore the SF ■ J term in the perturbation of the Maxwell stress (i.e. gradients 
of the background field tensor (oc J) gentler than those of the perturbations). 

The simple "linear poking" of the field tensor may now be written in a manifestly covariant manner 

F^8J V = -^{B-ikfe. (34) 

Naturally, evaluation of the elastic response of the field is straightforward in the rest frame of the fluid where one has 
F^SJv = ~(v Alt k^) 2 £ M . The only difference with the non-relativistic analog is that the Alfven speed is now weighted by the 
relativistic enthalpy of the fluid QV AU = ^F^F^. 

We are now all geared up to put together the pieces of the puzzle. In terms of the Lagrangian displacement vector field, 
the r.h.s.'s of Eqs [25 & 26] are to be balanced by the elastic response of the field tensor to the poking by £, c.f. Eq [34] 
(note that the radial scales of the transformation into the rest frame of the fluid cancel one another). This balance is locally 
equivalent to a M = — g^C, i.e. the covariant components of the fluid's acceleration respond to a force proportional to the 
displacement vector (with the unnormalized "spring constant" = {v AU k-^) provided by the field). By construction, both of 
these vectors are orthogonal to U and collinear. Furthermore, since £ M = (i'jcr) 2 ^ and 70- is a world scalar to be identified 
with the true co- moving frequency (as measured by an observer riding along with the fluid), it follows that £ v = («70") 2 £*. 
With these relations and the aforementioned equations for the tidal and coriolis terms, one arrives to lengthy component 
equations for £ r and ^ , for general f2 = U v /U t and negligible radial flow. 

In the case of circular geodesic flow, the equations simplify beautifully (horizontal regime) 
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^+ 2 7^Mr 3/2 ±a-2r 1/2 ) ? = (35) 

These immediately yield the sought after dispersion relation near a rotating hole 

( 7 a) 4 - [4 7 2 f^ (C± - \V) + 2 ? |] ( 7 a) 2 + 9 | [g| - 4 { § 7 2 D Q 2 ± }] = (36) 

where C± = 1 — § ± ^§72 corresponds to the C function of Novikov & Thorne (1973) for prograde orbits. 

Factoring out the extrinsic 5 dynamical frequency, fi± , one arrives to the normalized dispersion relation (with 7 u = fi± <r) 

* 4 -[C+X^ 2 +^C+*A]=8 (37) 
where 

A=-{§ 7 2 X>} and ±1 = 4 7 2 (C± - |X>) 

denote the normalized shear parameter and (co-moving) epicycle frequency (note that i x corresponds to the well known 
result of epicycle frequency as measured at asymptotic infinity). One thus sees that with the proper generalizations of the 
epicycle frequency and shear parameter, the local dispersion relation is identical with the Newtonian case in the limit of no 
fluid compression and lg • £ = (i.e., the "classical" Balbus-Hawley instability, Eq[ll]). 

Using the relation 7 2 = (1 ±a/r 3 ^ 2 ) 2 C^ 1 for cold, circular, geodesic flow (Novikov and Thorne 1973), one finds the fastest 
growing modes to conform with 

*--***--('±;j) 4 H5}' m 

which remains finite and close to the Newtonian value of i| for all radii outside the ISCO (and for any value of the rotation 
parameter) . 

To attach meaning to the polynomial functions that appear naturally in the dispersion relation for the magnetorotational 
instability, recall the range of radii that define particle dynamics in the Kerr geometry (Bardeen et al. 1972): 

(i) The marginally stable circular orbit (a.k.a. the ISCO), r ms , corresponds to the root of \± = 0. 

(ii) The radius of the circular photon orbit, r p h, is where C± = 0. 

(iii) The event horizon, r+, happens at the outer root of T> = 0. 

One therefore has the following ordering of radii for any value of the rotation parameter a: r ms > r p h > r+. As remarked 
by Bardeen et al. '72, when a = 1, the proper radial distance between these radii is non-zero in spite of "coinciding with the 
horizon", i.e. in spite of laying at the same Boyer-Lindquist radial coordinate. 

Inspection of Eq [38] now shows that — > + as r — > r+ h so the most unstable MRI modes go to large scale just outside 
the photon orbit. Moreover, utilizing that expression for q 6 in the unstable root of the dispersion relation, one finds the 
growth rate (or frequency!) to be given by 



„ 2 f 3 V } 2 [/ , a V 8C± / , „ a a 2 , a 3 a 4 \ /4C± 



(39) 



For a non-rotating hole, 

3 T> 

— > @ r = r p h(l + |), and the local growth rate a — —— > 2 

while for a rotating hole, the MRI quenching radii (for fastest growing modes) also occur just outside the circular photon 
orbit and may be readily extracted from the above relations. In Figs 1 and 2, we plot the general relativistic modifications the 
fastest growing linear wavemodes, wavenumbers and growth rates respectively, as functions of radius and for different values 
of spin parameter a. 

To go beyond this point, one would need to address global effects arising, for instance, from field curvature terms (see, 
e.g., Curry & Pudritz 1995, Ogilvie & Pringle 1996) and from the non-negligible radial velocity profile. Further investigation 
of the nature of the global instability is beyond the scope of this paper. 



4 DISCUSSION 

The MRI-in its most simple, local, incompressible variant-is found to operate virtually unabated down to the marginally stable 
orbit for massive particles. This radius is nearly coincident with the putative inner boundary of standard, thin accretion disks 
in the Kerr geometry. A vanishing epicycle frequency at r ms means that the fastest growing wavenumbers tend to be of a bit 

5 As defined, fi = U v /{/* reflects motion as observed in the Boyer-Lindquist frame, i.e., in a frame extrinsic to the fluid. It follows that 
the timescale associated with does not reflect a proper dynamical timescale. 
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smaller scale, g 2 : y| _ > 1 — Q(ar 3,/2 ), while growing faster than classically, id : § — » 1 + C(ar 3 ^ 2 ). The effects of strong 
gravity become truly significant only in a regime where circular, cold, geodesic flow is unstable (i.e., where x± < 0). 

Recall that particle trajectories with U v /U* = Q± exist inside r ma and all the way down to r p h but, in the presence of 
turbulent velocity fluctuations, body forces such as a radial pressure gradient would be required to confine the flow to such 
circular orbits. Although very little is concretely known about the accretion flow inside r ms , two rather robust remarks may 
be ascertained: The flow inside r m s cannot be supported centrifugally and it must therefore deviate from a standard thin disk. 
In addition, depending on the timescale for infall, the flow may not have time to cool significantly and advection of entropy 
will become progressively more important as r+ is approached. A robust prediction of this paper is the expectation that free 
energy tapping from the differential shear How goes on in the region immediately below r ms . 

One may envisage the situation inside the ISCO to evolve from a Mildly Advective Accretion Flow (MAAF) to a fully 
Advection Dominated Accretion Flow (ADAF) as the photon orbit is approached. In fully or partly advective accretion flows, 
such as those modeled by Popham & Gammie (f998), the angular velocity profile "peaks" precisely at r p h and quickly drops 
therein to match the angular velocity at r+. More importantly, when cooling by advection of entropy is moderately important- 
say, for advection fractions / ~ a few percent-the angular velocity profile departs very slowly from circular geodesic flow, 
U v /U* ~ f2± down to a region below the marginally bound orbit. The transition from nearly Keplerian to plunging orbits 
can be clearly seen in one of the very few global slim disk models where the cooling fraction is calculated explicitly: the 
fD models of Popham et al. (f998, albeit in the exotic scenario of a hyper-accreting black hole). In these models the radial 
velocity component is non- negligible when compared with the local sound speed (the sonic point generally occurs below r ms , 
even near r mb for low values of a), but v r is generally smaller than v v down to the region below r mb - (Note that the radial 
speed in the corotating frame, e.g. Gammie & Popham's (1998) V, is related to the same in the locally non-rotating frame 

The major limitation of the work presented herein is the presumption of negligible radial flow which greatly simplifies 
matters from the onset (see Eq [16]). At this point, it is unclear how much the results will change when full consideration 
is made for the radial inflow. Since the changes could be qualitatively significant-recent reports negate the reversal of the 
centrifugal force when the radial speed overwhelms the azimuthal component (Mukhopadhyay & Prasanna 2001, Prasanna 
200f)-this point should be the subject of close scrutiny in a future paper. Meanwhile, the adoption of an angular velocity 
profile corresponding to circular equatorial geodesic orbits seems a reasonable rough approximation in view of the above 
observations of advective flows. In this spirit, we argue below that the natural evolution of the MRI inside the marginally 
stable orbit is at least consistent with this assumption. 

Assume, in quasi-linear fashion, that the time- and length-scales provided by the linear dispersion relation reflect the 
growth and size of the dominant turbulent eddies to within factors of order unity to a few. Provided that o's/, simple linear 
growth/non-linear decay arguments (e.g., Araya-Gochez 1999) can be used to predict a predominantly toroidal field topology: 
The MRI constantly promotes radial/azimuthal field growth from "horizontal" velocity fluctuations, £ r ~ — while the 
coherent, background azimuthal shear flow converts this field into toroidal field at twice the rate of radial field generation. 
In the rest frame of the fluid, the tapping of free-energy associated with the shear flow becomes very rapid as the flow turns 
relativistic. Indeed, a co-moving observer measures the shear parameter, 2A/f2± to be §7 2 2? ~ § (1 ± a/r 3 ' 12 ) 2 V/C±, higher 
than the "Keplerian" frequency associated with the global dynamical timescale as seen at large distances (the redshift factor 
comes in because we chose to measure the angular frequency in terms of Boyer-Lindquist coordinates). 

The ratio T>/C± represents a gauge of the relative strength of two inertial terms, shear and coriolis. Setting aside the 
issue of radial flow for a moment, our dispersion relation suggests that as material approaches the region just outside of the 
photon orbit where C± vanishes, the slow branch of the dispersion relation (i.e. the MRI) is stabilized by the predominance 
of shear over the coriolis terms. Recall that the location of the circular photon orbit is the place where the centrifugal force 
reverses its direction: Inside r p h, increasing the velocity of a test particle pulls it in further (see, e.g., Abramowicz & Prasanna 
1990 and references therein). The limit of g fl — > + means that what was essentially a local instability becomes a global 
phenomenon. Although such a regime is formally outside the scope of the local analysis, one can anticipate a few rather 
interesting qualitative consequences. 

At first glance, the dispersion relation Eq [37] shows the appearance of an interchange, radially buoyant mode (T. Foglizzo 
priv. comm., Araya-Gochez 1999). More likely, this would simply imply the need for a steep radial stratification profile. Indeed, 
if the coherence lengthscale of the field were to reach the comoving length associated with the radial scale, Ar/VT), the disk 
could make a transition from centrifugally driven to magnetically driven: MRI modulated dynamics guarantee that the 
Alfven speed associated with the toroidal field at this large scale would be comparable to the orbital speed. Moreover, the 
field generated at large scales is less prompt to decay through reconnection and also more buoyant. This has very important 
consequences for the energy fraction going into-and persisting in-electromagnetic channels. 

The radial velocity profile will very likely change the expected outcome once the radial velocity becomes supersonic 
or super- Alfvenic, but some of the qualitative features of the this analysis may carry over when the full problem is solved, 
analytically or otherwise. If so, in this part of the so-called "plunging region" of the flow, the turbulent eddies will tend to grow 
larger while the field direction will tend to track the surfaces of null angular velocity gradients (no longer purely toroidal). 
The implied field topology is that of large-scale horizontal field domains. 
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Figure 1. Normalized wavenumber, q^, as a function of radius (in gravitational radii) for several values of the spin parameter a. 
Diamonds indicate the location of the marginally stable orbit, \ = 0, and triangles, the location of the marginally bound orbit. 

4.1 effects of radiation stress and neutrino trapping 

Precise assessment of the dynamical role of radiation in the general relativistic regime is hampered by the breakdown of one 
key assumption made to simplify the "linear poking" on the Faraday field tensor: Use of the enthalpy weighted specific four 
magnetic field in Eq [34] . On the other hand, one expects a photon gas-semi-contained by a neutral plasma through Compton 
scattering-to comprise a rather funny MHD fluid where the magnetic field is truly frozen only to the co-moving volume 
associated with the mass density but for which pressure perturbations do not behave adiabatically. It follows that when the 
fluid transitions into radiation pressure domination, compressive modes (e.g. toroidal field, non-axisymmetric modes) may lose 
pressure support in an unfavorable range of wavenumber phase-space (Agol & Krolik 1998). One can prove that the MRI falls 
squarely in such radiative heat conduction damping regime (Blaes & Socrates 2001) . Araya-Gochez & Vishniac (2002) show that 
the behavior of the energy equation is in some (algebraic) sense "quasi-adiabatic" for exponentially growing, non-propagating 
modes. Mathematically, this means that a real, analytical, slow-varying function of the scale of the perturbations, r(ik 2 /fe°), 
can be used to treat the energy equation in quasi-adiabatic fashion. Radiative heat conduction isotropizes the modes and, to 
zeroth order, one can use such quasi-adiabatic index in Eq [12] to anticipate that the effects of radiative heat conduction out of 
compressive toroidal modes is to increment the threshold of shear parameter where q & — > + from —2 to —A — > 1 + 2D/(1 + A). 
Nevertheless, since A cx T)/C± and C± — > @ r p h, the increase in shear threshold in this setting is rather inconsequential. 

Note further that the qualitative nature of energy deposition in radiation pressure dominated fluids is insensitive to the 
details of the (global) cooling but it is explicitly sensitive to the optical thickness of the relevant eddies. Thus, upon the onset 
of neutrino trapping in the neutrino cooling regime of hyper-accreting black holes, one may reasonably expect MRI modulated 
dynamics at p„ ^Pr + g (gas and radiation are tightly coupled) to resemble the standard disk case when p ra( j fep g as- Turner et 
al. 2001 report that the non-linear outcome of the MRI in this setting is a porous medium with drastic density contrasts as 
to cheat the Eddington limit at high accretion rates. Under nearly constant total pressure and temperature, the non-linear 
regime shows that density enhancements anti-correlate with azimuthal field domains (just as expected from the linear theory) 
and that turbulent eddies live for about a dynamical timescale while mass clumps are destroyed through collisions or by 
running through a localized region of shear. 

Since the turbulent eddies in the disk are largely instabilities of the toroidal field (at moderate values of the field), large- 
scale horizontal field domains near the marginally bound orbit would naturally force the baryonic component of the accretion 




© 2002 RAS, MNRAS 000, 000-000 



12 R. A. Araya-Gochez 



2 
1.8 
1.6 

d 1.4 
1.2 
1 

0.8 

1 2 3 4 5 6 7 8 9 10 
R [GM/c 2 ] 

Figure 2. Normalized growth rate, ct, as a function of radius for several values of the spin parameter a. 

flow into spatially segregated, massive clumps that occur near the nodes of non-axisymmetric (toroidal) MRI eddies (Araya- 
Gochez & Vishniac 2002) . This expectation motivates the picture of massive clumpy accretion suggested in the introduction. 



5 ENDING NOTES 

In summary, this work shows that the MRI is virtually unaffected by strong gravity outside the innermost stable circular 
orbit. Secondly, it indicates that the instability becomes non-local inside this region. Indeed, the MRI may leave behind a 
large scale, ordered field as the fluid heads in towards the circular photon orbit (with an orientation that tracks surfaces 
of null angular velocity gradients). Assuming incompressibility and the angular velocity profile of circular geodesic flow, the 
fastest growing modes die off while going to large scales at a radius just inside the marginally bound orbit. Accountability of 
compressibility as required to address the effects of radiation stress will bring the critical MRI quenching radius in, slightly 
closer to the photon orbit. Radiation stress, when significant, will diminish the growth rate while increasing the threshold of 
shear parameter to quench the MRI. 

Radial inflow will affect the global field topology but the details depend on poorly understood fluid trajectories in a 
region where cold, circular geodesic flow is unstable. As was pointed out by Krolik (1999), the standard assumption of 
ballistic orbits is never self-consistent for ideal MHD accretion inside r ma . Indeed, when magnetic turbulence is the culprit of 
angular momentum transport in the disk, the magnetic field energy density must become comparable to the rest-mass energy 
density of the fluid in the plunging region. Yet, unlike Krolik's suggestion, we do not believe that linear Alfven waves could 
efficiently transport energy from inside r ms ; the magnetic field there is still highly unstable, and the range of stability of such 
waves is limited by inertial forces. 

On the other hand, in this paper we demonstrate that energy deposition and angular momentum transport through the 
MRI go on virtually unscathed in the region just below r ms . An important note is the promptitude of this process at near 
Eddington rates since the MRI directly feeds the photon bath through compressive damping of the modes (Araya-Gochez & 
Vishniac 2002). Energy deposition into the radiation field thus occurs on the MRI timescale! On the other hand, near r p h 
the flow will inevitably transition into advective cooling. Assessing the magnetic field dynamics in the region v ms ^> v ^ 7*mb 
is essential to predict the efficiency of accretion, and to address some large scale effects such as jet launching and disk-hole 
coupling. 

At highly super-Eddington accretion rates (such as those expected in the prompt stages of hyper-accreting black hole 
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formation), the fluid may possess non-trivial amounts of internal energy per unit rest mass of baryons. For such a hot 
MHD fluid, r ms does not represent a significant boundary to the disk/flow and such may occur rather closer to r p h. This 
stresses the importance of addressing MHD processes in the region above the circular photon orbit. Along these lines, we 
have motivated the provocative conjecture that copious gravitational wave losses ensue through black hole ringing when a 
hyper-accreting black hole enters the accretion regime where neutrino trapping occurs. This argument, which combines linear 
regime phenomenology with the latest numerical results from accretion in radiation-stress dominated environs, leads to a 
picture of near-hole accretion where large-scale horizontal field domains channel the flow into massive clumps that "thump" 
the hole. 

Lastly, note that a strong, toroidal field topology is ripe ground for MHD instabilities that promote poloidal field generation 
such as the Parker and radial interchange instabilities (in the vertical and horizontal regime respectively). This instabilities 
could provide a physical justification for desirable field topologies invoked in jet launching and the Blandford-Znajek processes. 
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